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We investigate a new hybrid-model universe containing two types of dark 
matter, one "warm" and the other "hot". The hot component is an ordi- 
nary Ught neutrino with mass ~ 25/1^ eV while the warm component is a 
sterile neutrino with mass ~ 700/i^ eV. The two types of dark matter arise 
entirely within the neutrino sector and do not require separate physical ori- 
gins. We calculate the linear transfer functions for a representative sample 
of warm-plus-hot models. The transfer functions, and results from several 
observational tests of structure formation, are compared with those for the 
cold-plus-hot models that have been studied extensively in the literature. On 



the basis of these tests, we conclude that warm-plus-hot dark matter is essen- 
tially indistinguishable from cold-plus-hot dark matter, and therefore provides 
a viable scenario for large scale structure. We demonstrate that a neutrino 
mass matrix can be constructed which provides the requisite dark matter 

constituents, while remaining consistent with all cosmological bounds. 
PACS number(s): 12.60-i, 13.15.+g, 14.60.Lm, .Pq, .St, 14.80.Mz, 98.90.-k, 

.Cq, 95.35.-Fd 
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I. INTRODUCTION 



Growing evidence against both the cold dark matter (CDM) and hot dark matter (HDM) 
cosmological models has prompted researchers to turn to alternative scenarios. Mixed dark 
matter (MDM) cosmologies [1,2], in which the universe contains an admixture of HDM, 
CDM, and ordinary matter, have captured most of the attention. By properly tuning the 
amount of HDM in the mix (20% — 30% in most scenarios) one can find models with a 
balance of small {5h~^ Mpc) and large (25 — 50/i~^ Mpc) scale power consistent with current 
observations [2] {h is the Hubble constant today in units of 100 kmsec~^Mpc~^). Recent 
reports [3] from the Los Alamos National Laboratory suggesting that i^e and/or i/^ have 
mass in the few eV range (just what the MDM models demand) have heightened interest in 
these models. 

The essential property of HDM is that relativistic dark matter particles free-stream out 
of high-density regions, and therefore fluctuations are damped on scales smaller than the 
horizon size at the epoch when the HDM becomes non-relativistic. For ordinary neutrinos, 
the canonical HDM candidate, the size of the first objects to form is Aps = 13Mpc (fi^h^)"^, 
where fl^h^ = mu/QSeV is the energy density in dark matter divided by the critical density, 
and nil, is the neutrino mass. In MDM models, the cold component clumps on all scales. 
By shifting some of the mass density from hot to cold matter, one boosts small-scale power 
relative to large-scale power. This is the essence of the MDM models. 

The greatest weakness of MDM models is that they require two types of dark matter 
with comparable mass densities. The hot component is usually taken to be an ordinary 
neutrino while the cold component is thought to be one of the standard CDM candidates 
(e.g., the hghtest SUSY particle or the axion). A single dark matter candidate requires "new 
physics" beyond the Standard Model of strong and electroweak interactions. At first glance, 
a second type of dark matter would require adding another new sector to the particle physics 
theory, with the unexplained coincidence that the two types of dark matter give comparable 
contributions to the total mass density of the universe. (See, however, ref [1] for an attempt 
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to provide a particle physics connection between hot neutrinos and cold axions.) 

In this paper, we explore the possibility that both the hot and cold dark matter compo- 
nents of a MDM universe are neutrinos. We propose that one of the active (weak interactions 
with matter) neutrinos is the hot component of the dark matter while one of the sterile neu- 
trinos (no weak interactions) is a warm component. For the sterile neutrinos not to be hot, 
their mass to temperature ratio, ms/T^, must be greater than that of ordinary HDM. We 
therefore require that the sterile neutrinos decouple from the rest of the plasma prior to the 
electroweak (EW) phase transition. During the subsequent EW and QCD phase transitions, 
the sterile neutrinos cool relative to the coupled plasma (including the active neutrinos). 
is lowered and, assuming a fixed energy density, is higher. ms/Tg will be roughly 15 times 
higher than that of a standard (~ 25/i^ eV) HDM neutrino. This type of particle is usually 
referred to as warm dark matter (WDM) and was first considered in the early 1980's [4]. 
Large scale structure in a pure WDM universe with m/T this large is much like large scale 
structure in a CDM-dominated universe. Even if the sterile neutrinos decouple between the 
QCD and EW phase transitions (m^/Tg roughly 8 times that of standard HDM neutrinos) 
they still behave much like cold dark matter [5]. 

Our main motivation is to show that two types of dark matter can arise simultaneously 
from the neutrino sector and provide a viable scenario for structure formation. We consider 
it a positive feature of the scenario that the two types of dark matter have masses which differ 
by only an order of magnitude or two. However, readers familiar with the standard lore of 
neutrino mass generation, will recognize our neutrino mass spectrum as unconventional. In 
the usual see-saw mechanism for generating small masses for active neutrinos, right-handed 
fields are typically very heavy due to the presence of large Majorana mass terms Mr (~ 
lO'^^^^GeV). The Dirac mass terms, which couple left and right-handed fields, are assumed 
to be comparable to the masses of the associated charged leptons or up-type quarks, i.e. 
mz) ~ MeV-GeV. Assuming no left-handed Majorana mass terms, the physical eigenstates 
corresponding to the active neutrinos have masses ~ (m|,/A'//j), provided rriD <^ Mr. While 
the see-saw mechanism may be the most attractive scenario for explaining small (but non- 
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zero) masses for the active neutrinos, there is no a priori reason for all three right-handed 
fields to be so heavy. The models constructed later are counterexamples where the physical 
state corresponding to one of the sterile neutrino types {i.e. mostly right-handed fields) has 
mass rus ~ 700/i^ eV. 

II. LARGE SCALE STRUCTURE WITH TWO TYPES OF NEUTRINO 

Our model universe consists of four components: two types of massive neutrinos, two 
massless neutrino species that are treated as a single component, and a mixture of photons, 
baryons, and electrons that is treated as a single component ideal fluid. We assume primor- 
dial perturbations that are adiabatic, Gaussian and scale-free with a spectrum Pp{k) = Bk"', 
where k is the wavenumber of the perturbations measured in units of Mpc~^. The COBE 
DMR experiment probes energy density perturbations on very large scales where there is 
little modification of the primordial spectrum. It finds B = 8.2 x 10^/i~*^Mpc^ and estimates 
of n that are consistent with n = 1 [6,7]. We use this value of B where normahzation of the 
power spectrum is required and set n — 1. We also set h — 0.5 and ignore baryons. h — 0.5 
is the value used in most studies of CDM and MDM [2,8] though it may be in conflict with 
measurements of the Hubble constant [9] . We emphasize that our main purpose here is to 
compare warm-plus-hot models with currently popular models such as MDM and not to 
carry out a detailed comparison with the observations. 

Since the perturbations at early times are small, one can follow the initial stages of their 
evolution using linear theory. In fact, linear theory is adequate for studying structures on 
all scales significantly larger than Sh"^ Mpc. Quantitative tests on smaller scales probe 
non-linear structures, and therefore require N-body simulations. 

The active neutrinos decouple from the photon-baryon plasma when they are still rel- 
ativistic and when the photon temperature Xy ~ 1 MeV. Their background distribution 
function is therefore /o(p) = {eP/^-^ + where = (4/11)^/^7^, p ^ {El - mlf'^, and 
El, is the neutrino energy. We assume that the background distribution function for the 
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sterile neutrino species also has a Fermi-Dirac shape; fo{p) — (3 [eP/'^'' + ij . This distri- 
bution function is appropriate for a neutrino that decouples when it is relativistic provided 
(3—1 and Tg — (4/11)^^^ T^, where is the effective number of relativistic de- 

grees of freedom in the universe at the time of decoupling. The masses for the two neutrinos 
satisfy the relation 

+ (p) = 93n,/i2gv. (1) 

For a particle that decouples between the EW and QCD phase transitions, ~ 61, whereas 
a particle that decouples prior to the EW phase transition has g^ ~ 107. We refer to models 
of these types as WPHl and WPH2, respectively. 

The above discussion assumes that the sterile neutrinos were coupled to the baryon- 
photon plasma at some early epoch. One possibility is that coupling occurs at very high 
energies where unknown physics (e.g. from Grand Unified Theories) operates. Alternatively, 
oscillations between active and sterile neutrinos may be responsible for bringing the sterile 
neutrinos into equilibrium [10,11]. In addition, Dodelson and Widrow [12] have shown that 
under certain assumptions, the momentum space distribution function for sterile neutrinos 

1 /3 

produced through oscillations has a Fermi-Dirac shape with Tg = (4/11) T^, but /5 < 1. 
These models lead to exactly the same large scale structure phenomenology as the early- 
decoupled particle models [5] and so our discussion of structure formation applies to both. 

It is convenient to define the linear transfer function T{k) = {6p{k)/p) / {Sp{k oo)/p) 
which describes the evolution of the power spectrum through the linear regime. The evolved 
linear power spectrum is then P{k) — Bk^T{kY. The hnear transfer functions for a variety 
of models are shown in Figure 1. As we can see, the WPH2 models are very close to the 
MDM models; the WPHl models have a more distinct shape, though even in this case the 
differences are likely to be observationally indistinguishable. A more quantitative comparison 
of the models is made by calculating various integrals (with appropriate window functions) 
of the power spectrum. In particular, we calculate erg, EP (the excess power on 25 /i~^Mpc), 
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and (7o.5. Here, ctl is the mass excess on Lh 



Mpc defined by the integral 



1/2 



(2) 



where W{x) — 3(sinx — xcosx) /x^ is the top hat window function. Fluctuations in the 
mass density are related to fluctuations in optically-selected galaxies through biasing, where 
ho(TL gives the fluctuation in these galaxies on the scale L/i~^Mpc. hoi the optical biasing 
parameter, can depend on scale, though in the simplest models it is assumed to be constant. 
Davis and Peebles [13] find that hoOs ^ 1 and therefore l/us is a measure of the optical bias. 

As mentioned above, CDM falls short because it has too little power on large scales 
relative to small scales. Wright et al. [6] introduce the quantity EP = 3.4(725/(78. This 
definition is such that EP — 1 for standard CDM, whereas consistency with the APM 
angular correlation function requires EP = 1.30 ± 0.15. Finally, cro.5 gives a rough (and 
probably overly pessimistic) approximation for 1 + Zg{, where Zgi is the redshift for galaxy 
formation. (Alternative methods for estimating can be found in ref. [14] and references 
therein.) 

The results for ag, EP, and (70.5 for a variety of models are given in Table 1. For MDM, 
CDM, and HDM, we use analytic fitting functions for T(k) taken from Holtzman [15]. We 
see that the WPHl, WPH2, and MDM models with a 30% hot component give similar 
results for EP and (Tg. It has been argued [2] that MDM 30% satisfies large scale structure 
tests on 8 — 25/i~^Mpc scales, and therefore these conclusions should hold for our models 
as well. We do note, however, the results for (Tq.s in Table 1, and also discussion in the 
literature (e.g. ref [2]), suggest that any model with a significant hot component may have 
trouble if observations push the redshift of galaxy formation Zgf > 5. This situation can be 
improved by decreasing the amount of hot matter in the mix, though this also leads to a 
decrease in power on large {25h~^ Mpc) scales. Since galaxy formation depends on nonlinear 
and nongravitational physics (e.g., hydrodynamics) the above conclusions should be used 
with caution. 
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Model 






0-8 


EP 


0-0.5 


HDM 


23 


— 


0.87 


1.46 


1.0 


WPHl 20% 


4.6 


100 


1.00 


1.18 


2.1 


WPHl 30% 


7.0 


90 


0.93 


1.27 


1.6 


WPH2 20% 


4.6 


190 


1.02 


1.17 


2.6 


WPH2 30% 


7.0 


160 


0.94 


1.26 


2.0 


MDM 30% 


7.0 


> 10^ 


0.96 


1.30 


2.6 


CDM 




> 103 


1.48 


0.94 


9.6 



Table 1 Results for ag, EP, and ctq.s for a variety of models. The results are based on linear 
theory calculations. Percentages in column 1 for the hybrid models indicate fraction of dark 
matter in the mix. 

III. NEUTRINO MASSES - MODELS AND CONSTRAINTS 

At temperatures comparable to Tqcd, sterile neutrinos interact with ordinary matter 
by mixing with active neutrinos. From detailed studies of neutrino oscillations in the early 
universe [10,11] the conditions that the sterile neutrinos are out of equilibrium below Tqcd 
(and therefore do not spoil the predictions of Standard Big Bang Nucleosynthesis [16]) is 

(5m2 < 3.6 X 10~^(sin2 2^o)-2 eV^ ; for sin ^ 10"^ , (3) 

where dm? = m'j.—m'^ is the mass squared difference between the states, and 9o is the vacuum 
mixing angle between them. This is provided that the energy of the states E ^ V 6m?, so 
that they do not decohere. However, for WPHl it does not matter if oscillations bring 
the sterile neutrinos into equilibrium above Tqcd- As long as the sterile neutrinos decouple 
prior to the QCD phase transition, entropy transferred from the quark-gluon plasma into the 
interacting gas will dilute the sterile neutrino number density relative to the active neutrino 
number density. For this reason there is no constraint for sin ^ 10"^ (The exact limit 
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depends on the adopted value for Tqcd)- For WPH2, the demand that the sterile neutrinos 
are out of equihbrium below the EW phase transition at Tew ~ 300 GeV, is essentially that of 
Eq. (3) except that the region where no constraint applies becomes sin 6*0 ^ 10~®. (However, 
these bounds are model dependent. In the singlet majoron model [17], for example, they are 
greatly weakened due to the contribution of the majoron background to the active-neutrino 
self-energy [18].) 

A specific mass matrix which is useful for further discussion is 
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M' 



































(4) 



where we label our weak eigenstates as 



(L) VfiiL) ^t{L) l^a{R) J^(3{R) ^■y{R) 



(5) 



There is no reason to identify a, P and 7, with e, and r. Mg^ is kept non-zero only to 
simplify later calculations. In (4) the matrix elements rriij and M^j correspond to Dirac and 
Majorana mass terms, respectively. As in most models, the upper left 3x3 submatrix of (4) 
is taken to be zero to be consistent with measurements of the width of the Zq (see however 
[19]). The other vanishing terms of (4) are set to zero for clarity, and in general they need 
not be identically zero for our scenario to remain viable (see later discussion). 

Diagonalization of (4) leads to six mass eigenstates, which in the limit rriij ~ m <S M^j ~ 
M have mass eigenvalues 

,2\ 



m 



mi = 0; 7712, tTis = O \ -j-^ 1 ; 777.4 = M'; 7775, 7776 — ^ 



(6) 



For judicious choices of the entries in (4), the physical state corresponding to either 1712 or 
7773 will correspond to the (active) HDM component and that corresponding to 7774 will be 
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the (sterile) WDM component. Those corresponding to m.5 and me can be made sufficiently 
heavy with decay rates (into some combination of light neutrinos and possibly hght scalars) 
short compared to the age of the Universe. The mass eigenstates corresponding to 1712 and 
ms will be a mixture of the t'e(L), and i^t{l) states but with a small admixture 0{m/M) 
of the sterile and v^(r)- lim/M were large enough, we see from (4) that the constraint 
w? /M < 5 X 10~^ eV would be imposed. Since w? /M is the mass of our light active state, 
this would imply that our HDM candidate could not possess a cosmologically interesting 
mass. However, if we take M ~ 1 TeV, then E < ^J{5w?). As such, no oscillations will 
occur and we are relieved of this mass limit. The mixing angle with the lighter sterile state 
i>'4 is zero so there is no oscillations into this state. Prom this we see that a mass matrix 
with m ~ 3 MeV, M' ~ 200 eV and M ~ 1 TeV, would satisfy the oscillation constraints 
and yield the desired neutrino mass spectrum. 

The mass matrix (4) has been chosen to elucidate our discussion. The first key feature 
of our scenario is the introduction of a new mass scale M' <^ M. M may be the vacuum 
expectation value of one Higgs field and M' that of another. Alternatively, there may be 
only one new Higgs with different Yukawa couplings to the right-handed neutrino states. 
The smallness oi M'/M is the main "unnatural" aspect of our scenario. However, to put it 
into perspective, the above-quoted value for M'/M ~ 10~^ is not enormously smaller than 
some of Yukawa couphngs in the Standard Model, such as g'y(e) = rUe/Myjeak — 10~^, and is 
much larger than the ratio of Higgs vacuum expectation values postulated in Grand Unified 
Theories, Mw/Mqut ^ 10"^^. 

A second feature of our scenario is the non-participation (or at least reduced participa- 
tion) of the light sterile state in the usual mixings of the left and right-handed neutrinos 
and in the consequent see-saw mechanism. In (4) this is accomphshed by the introduction 
of matrix terms identically equal to zero. These zeros may be enforced by the imposition 
of global U(l) symmetries. For example, in the singlet majoron model [17] the Dirac mass 
terms arise from the coupUng of the neutrinos to the the usual Standard Model Higgs dou- 
blet, H. In addition, one introduces a complex scalar (Higgs) field $ which is a singlet under 
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SU(2)xU(l) gauge transformations, but carries non-zero Lepton-number. $ gets a vacuum 
expectation value generally taken to be ($) GeV — lOTeV. The right-handed Majorana 
masses come from their coupling to ($). If one extends this model by adding two new scalars 
$1 and $2 (or by adding another sterile neutrino), one can easily enforce the explicit zeros 
of (4) without generating any others - for example using the Lepton number assignment: 

^'e(L) ^'/i(L) 1^t{L) ^a(K) ^I3{R) l^y{R) ^1 $2 H 

L 1 1 1 2 1 1 -2 -4 

In general the zeros of (4) need not be identically zero. That they be small relative to the 
other scales in the matrix would, in most cases, suffice. (The exact value required depends 
on the details of the model.) If indeed they are non-zero, then neutrino oscillations could 
be the mechanism that brings the sterile neutrinos into equilibrium at early times. There 
are two constraints we must impose on the non-zero mass terms. First, the admixture of 
active states in the heavier sterile mass eigenstate must be small enough that oscillations 
will not keep the sterile state in equilibrium below either Tew (for WPH2) or Tqcd (for 
WPHl). (An explicit calculation in the context of the singlet majoron model shows that U4 
is decoupled prior to Tew-) Second, experimental limits on neutrino-less /3/3-decay require 
that [20] 

{mj = |S,C/>.J2 < 1.1 eV, 

where Uij is the unitary matrix diagonalizing (4). This imposes limits on the values of rriij. 

The discussion above provides an existence proof of our scenario. The mass matrix 
introduced and the constraints imposed on it are independent of the underlying particle 
physics model. However, there are other conditions we must satisfy, and to investigate them 
in any meaningful way requires the specification of a particular particle physics model. The 
model-dependent constraints are: 

(i) The decay of the light sterile state Us must be inhibited. The most general neutrino decay 
scheme can be written i/^ — > i/^ -|- X where i/^ is some lighter neutrino (in our scenario one 
of the active states), and X is some set of bosons and/or pairs of lighter fermions. 

(ii) Any new light degrees of freedom (e.g., the majoron) must be decoupled prior to Tqcd, 
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so as not to contribute significantly to the energy density of the universe at the epoch of 
nucleosynthesis. To ensure this, we assume that our new particle does not couple in any 
significant way with the active neutrino states or with ordinary matter. If it couples strongly 
only with the sterile neutrinos, then it to will be diluted away by entropy transfer at the 
electroweak and QCD phase transitions. Again, in the context of the singlet majoron model, 
this proves to be the case for the values of m, M' and M quoted. 

(iii) The heavy degrees of freedom (here v^^) must annihilate or decay away before they 
can dominate the energy density of the universe. Again, in the singlet majoron model, the 
heavy states decay away on a sufficiently short time scale. 

(iv) new /3/3-decay modes must be adequately suppressed. In the singlet majoron model, 
Qeej the effective couphng of the Majoron to must be < 7 x 10~^ [21]. This is satisfied by 
the model presented here. 

The above constraints are clearly model dependent in that they rely on the details of the 
theory which ascribes mass to the neutrino sector. Any detailed investigation of the validity 
of our scenario therefore requires knowledge of the underlying particle physics. We have 
sometimes highlighted the possibilities for our scenario in the context of the singlet majoron 
model. Many other models possessing sterile neutrino fields, including more complicated 
majoron models, can also be constructed. 

IV. CONCLUSION 

We have investigated the possibility that the universe possesses a mixture of warm and 
hot dark matter. The hot component is identified with a known light neutrino with mass ~ 
25/i^ eV; and the warm component with a sterile neutrino with mass ~ 700/i^ eV. These mass 
ranges satisfy the normal cosmological mass limits for stable neutrinos, provided the sterile 
neutrinos decouple at an early epoch. From calculations of the linear transfer functions for 
such a hybrid model, we show that mixed warm-plus-hot dark matter is a viable cosmological 
scenario, and in most respects provides a better fit to observational data than either the 
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standard HDM or CDM models. It is also interesting to note that the 700/i^ eV sterile 
neutrino can evade the usual phase space constraints [22] . 

The introduction of our warm dark matter particle has come at a price, namely, the ne- 
glect of theoretical prejudice that all right-handed neutrino fields be very massive (m >GeV). 
We emphasize, however, that a priori the mass terms of matrix (4) are unknown, and there 
is no fundamental reason that sterile neutrinos cannot possess masses in the eV range. 

In addition, we have discussed how such particles satisfy all standard cosmological bounds 
given certain restrictions on their oscillation properties. If theoretical prejudice is relaxed 
just a little, we see that a solution to the problems of the formation of large scale structure 
could reside purely in the neutrino sector. 

We find it suggestive that one can accommodate a hybrid dark matter universe entirely 
within one sector avoiding the need for separate physical origins of the dark matter con- 
stituents. Given these arguments, and given the recent tentative experimental evidence for 
a hot neutrino component, we believe this model deserves further consideration. 
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